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ABSTRACT 

We report the results of a search for gas phase atomic metals in the circumstellar envelope of the asymptotic giant branch 
carbon star IRC+10216. The search was made using high resolution (A/AA = 50 000) optical absorption spectroscopy 
of a background star that probes the envelope on a line of sight 35" from the center. The metal species that we detect 
in the envelope include Nal, Kl, Cal, Call, Crl, and Fel, with upper limits for All, Mm, Til, Tin, and Srn. The 
observations are used to determine the metal abundances in the gas phase and the condensation onto grains. The metal 
depletions in the envelope range from a factor of 5 for Na to 300 for Ca, with some similarity to the depletion pattern 
in interstellar clouds. Our results directly constrain the condensation efficiency of metals in a carbon-rich circumstellar 
envelope and the mix of solid and gas phase metals returned by the star to the interstellar medium. The abundances of 
the uncondensed metal atoms that we observe are typically larger than the abundances of the metal-bearing molecules 
detected in the envelope. The metal atoms are therefore the major metal species in the gas phase and likely play a key 
role in the metal chemistry. 
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1. Introduction 

Metals are important constituents of the interstellar 
medium. Abundant metals such as calcium and iron are 
among the most heavily depleted elements in dense clouds, 
and therefore form a significant component of interstellar 
dust. The process of dust formation is not well understood 
(e.g., Draine 2009), but much of the raw material is supplied 
by mass loss from evolved stars, especially the asymptotic 
giant branch (AGB) stars. 

In the winds of oxygen-rich AGB stars, the dust is pri- 
marily in the form of metal silicates, so that the metals play 
an important role in the formation of dust and the mass 
loss process. In carbon-rich AGB stars, the situation is less 
clear. The dust is believed to consist mainly of graphite or 
amorphous carbon, and silicon carbide. The extent to which 
metals contribute to dust formation, and the form in which 
they are returned to the interstellar medium is essentially 
unknown. 

In this paper, we report the first comprehensive search 
for gas phase atomic metals in a carbon-rich circumstel- 
lar envelope. IRC+10216 (CW Leo) is the nearest carbon 
star with a thick circumstellar envelope, and serves as an 
archetype for the study of mass loss on the AGB. The star 
is relatively faint, ^16.0 mag. in the i?-band and much 
fainter at shorter wavelengths because of obscuration by 
the envelope, but at longer wavelengths the circumstellar 
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dust and gas are seen in emission, and are brighter than 
for any similar object. IRC+10216 has therefore been in- 
tensively observed, with more than 50 molecular species 
detected in the envelope, including several metal bearing 
species (e.g., Olofsson 2005; Ziurys 2006a). 

The technique that we use here to search for atomic 
metals in the envelope of IRC+10216 is optical absorp- 
tion spectroscopy, using a background source of illumi- 
nation. The star itself is too faint and its spectrum too 
complex to serve as a useful source for detailed study, al- 
though circumstellar C2 and CN have been observed in 
this way (Bakker et al. I1997p . There are, however, other 
stars in the field. IRC+10216 is nearby, at a distance of 
r~> 120 pc (e.g., Ramstedt et al. 2008), so the circumstel- 
lar envelope extends a considerable angle on the sky. It 
is detected out to a distance of ~ 3' from the center in 
millimeter CO emission (Huggins et al. 1988), and 9' in 
infrared dust emission observed with IRAS (Young et al. 
1993). Although IRC+10216 is at a relatively high galactic 
latitude (I = 221°, b = +47°), there are several stars in 
this region of the sky that are candidates for background 
sources, as seen in the wide field image in Fig. 1 of Mauron 
& Huggins (1999). 

One of these stars, Star 6 in the UBV photometric se- 
quence of the field by Mauron et al. (2003), is well suited 
for absorption line studies. It is located behind the enve- 
lope at an angular offset of 35" from the center, and is 
bright enough for high resolution spectroscopy. This star 
has been observed with the UVES spectrograph at the VLT 
by Kendall et al. (2002). Their main objective was to search 
for diffuse bands that might originate in the circumstellar 
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gas. No diffuse bands were found, but these authors noted 
deep absorption lines of Na I and K i, which they attributed 
to circumstellar gas. 

Here we report a comprehensive search for metal lines 
in the circumstellar envelope of IRC+10216 along this line 
of sight. Our objectives are to measure the degree of metal 
depletion in this carbon rich-environment, and to deter- 
mine the distribution of solid and gas phase metals returned 
by the star to the interstellar medium. Sect. 2 describes 
the observational material. Sect. 3 presents the absorption 
lines, and Sect. 4 the derived column densities and abun- 
dances. The results are discussed in terms of dust formation 
in Sect. 5, and metal chemistry in Sect. 6. Our main con- 
clusions are given in Sect. 7. 



2. Observations 

The observations of the envelope of IRC+10216 were made 
using absorption line spectroscopy with Star 6 (USNO 0975- 
0633-6975) as the background source of illumination. Star 6 
lies 35" from IRC+10216 at position angle 165°. Its visual 
magnitude is F=16.0, and it is the nearest star to the cen- 
ter of the envelope suitable for high resolution spectroscopy. 
The field is shown in Fig. 1, where the envelope is seen in 
dust-scattered Galactic light. From its magnitude and spec- 
tral type (type G), star 6 is at a distance of ~ 1400 pc, well 
beyond IRC+10216, as discussed by Kendall et al. ([2UtJ2")) . 

The observations were made with the UVES spectro- 
graph at the VLT, and were previously used by Kendall 
et al. (2002) to search for diffuse bands. The data were 
obtained over seven nights in December 2000 and January 
2001, with an effective exposure of 4 hr at each wavelength. 
The sky transparency was moderate to good, and the see- 
ing was 0'.'8. The data consist of 4 spectra, each covering 
the complete range from 3 000 to 10 000 A with a resolv- 
ing power of 50 000 (Av — 6 kms -1 ). These spectra have 
been reduced with the ESO UVES pipeline and have been 
summed after correction to the heliocentric reference frame. 
All radial velocities given in this paper are heliocentric, ex- 
cept where specified otherwise. The signal-to-noise ratio of 
the final spectrum at 4200 A is ~ 60 per resolution element. 

The heliocentric velocity of Star 6 measured from nu- 
merous photospheric absorption lines is +52.4+0.6 kms . 
For comparison, the systemic heliocentric radial velocity of 
IRC+10216 is —19.3 kms -1 , which is derived from millime- 
ter observations of the envelope (ulsr = — 26 kms" 1 , Loup 
et al. 1993). Thus a given line in the envelope of IRC+10216 
is well separated in velocity from the same line formed in 
the photosphere of Star 6. In addition, the expansion ve- 
locity of the envelope is 14.1 kms -1 with a small turbulent 
component (Huggins & Healy 1986), so the absorption lines 
in the envelope are expected to be wide (~ 30 kms -1 ). 
In contrast, the weak stellar lines are relatively narrow, 
with widths comparable to the instrumental profile, and 
the strong stellar lines have characteristic damping profiles 
with broad wings. 

Table 1 lists the metal lines that we searched for in 
the envelope of IRC+10216. Column 2 gives the laboratory 
wavelength in air (A ) for each transition, and column 3 
gives the oscillator strength (/) from Morton (1991, 2000). 
The line list includes the strongest ground-state transitions 
for the species that are potentially observable in the wave- 
length range covered. These lines are seen in interstellar 




Fig. 1. V^-band image of the circumstellar envelope of 
IRC+10216, made with the VLT. The field size is 90" x 90". 
The mass-losing carbon star is located at the center of the 
image, and Star 6 is the bright source near the bottom, 35" 
from the center. North is to the top, East to the left. 

clouds and/or in circumstellar environments with similar 
physical conditions. 

At long wavelengths, the spectrum of the background 
source, Star 6, is relatively free of photospheric lines, so 
an absorption line arising in the intervening circumstellar 
envelope is straightforward to identify. At shorter wave- 
lengths, the spectrum is more crowded with photospheric 
lines, and absorption by the circumstellar envelope is often 
more difficult to identify. In these cases we use a template 
technique to extract the envelope signal. We first make a 
least squares fit to the spectrum of Star 6 using a template 
spectrum covering a region of ~ 15 A around (but exclud- 
ing) the expected envelope line. We then use this as the 
effective continuum to search for residual absorption from 
the envelope. 

One form of template that we tried was based on the 
library of synthetic spectra from Coelho et al. (2005), but 
these did not produce a good match to the photospheric 
spectrum. A second approach that was successful, was us- 
ing a scaled solar spectrum. By chance the spectrum of 
Star 6 is very similar to that of the Sun, and a template 
based on scaling the solar spectrum gives a close match to 
the stellar spectrum. The solar spectrum that we use here 
was published by Dclbouille et al. (1973) and covers the 
wavelength region 3 000-10000 A with a step size of about 
0.0125 A. In order to fit the stellar spectrum -F*(A), we shift 
the solar spectrum Fq(X) to match the radial velocity, re- 
bin it to match the resolution, and scale the flux according 
to: 

F.(\) = aF e (\) + , 

where a and /3 are constants that are determined from the 
least-squares fit to the local region of spectrum under con- 
sideration. 

The use of this template as the effective continuum re- 
veals spectral features of the circumstellar envelope that are 
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Upper limits are 3(7 . 

D-lines are saturated, N based on the 3302 A line. 
Based on the 7699, 7665 A doublet. 
: indicates 2-3(7 features; N based on weighted mean. 
Based on the 3860, 3720 A lines which have better S/N. 



not otherwise easily observed. Examples of the application 
of this technique are given in Sects. 3.3 and 3.4. 



3. Envelope absorption lines 

3.1. Ki 

Fig. 2 shows the observed spectral regions covering the K I 
doublet. The 7665 A line is in the upper panel and the 
7699 A line is in the lower panel. Strong photospheric K I 
lines of Star 6 are present at the stellar radial velocity of 
V r = +52.4 kms -1 , which is marked with an arrow at the 
right of each panel. Also present are two weaker photo- 
spheric lines (marked +) and four telluric lines (marked 

The dominant feature in the middle of each spec- 
tral region is strong Ki absorption from the envelope of 
IRC+10216. This can be unambiguously ascribed to the 
circumstellar envelope. It is centered near the systemic ra- 
dial velocity of IRC+10216 at -19.3 kms" 1 (marked with 
an arrow in the figure) and has the broad width expected 
from the expanding envelope. With the low interstellar ab- 
sorption in this region of the sky [Eb-v & 0.03), any inter- 
stellar Ki would contribute a few percent of the strong 
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Fig. 2. Spectra of the Ki doublet at 7664.91 A (upper 
panel) and 7698.97 A (lower panel). In each panel, the 
arrow on the right marks the photospheric line (V T = 
+52.4 kms -1 ), and the arrow in the center marks the 
envelope absorption of IRC+10216 (V T = -19.3 kms -1 ). 
Other photospheric lines are marked +, and telluric lines 
are marked o. 



lines observed, based on the survey of interstellar Ki by 
Chaffee & White (1982). 

We also searched for the much weaker ground state K I 
line at 4044 A which has an /-value ~ 110 times less than 
the 7664 A line. The 4044 A line was not detected, and the 
upper limit is given in Table 1. 

Although the K I doublet lines formed in the envelope 
are broad, they are clearly composed of several components. 
The presence of the components is important because it af- 
fects the saturation of the lines (see Sect. 4). Assuming 
gaussian broadening, we found that a best fit synthesis 
taking into account the instrumental profile requires four 
components to match the line shapes. A preliminary fit 
with unconstrained parameters gave consistent estimates 
for the radial velocities of the individual components in each 
line. We then fixed the velocities at their mean values, and 
fit the profiles by varying the relative column densities of 
the components and their line broadening b-values, which 
we constrained to be the same for both lines. The final 
mean parameters are given in Table 2, where the strength 
is the column density of each component, relative to the 
strongest component. The uncertainty in the velocities is 
^ 0.5 kms -1 and the relative strengths found from the two 
lines are the same to within ^ 5%. Thus the results from 
the two lines are in good agreement, as expected from their 
similar line profiles. 
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Fig. 3. Comparison of the observed (black) and synthesized 
(grey) profiles of the K I 7698.97 A line, using the param- 
eters given in Table 2. The individual components are also 
shown. 



To illustrate the quality of the multi-component fit, 
Fig. 3 shows the results of a synthesis of the 7699 A line 
using the parameters of Table 2. It can be seen that the 
model provides an excellent fit to the observational data. 

We interpret the multiple component profiles of the K I 
lines as the effects of the multiple shell structure in the 
envelope (Mauron & Huggins 1999, 2000; see also Fig. 1), 
where the line of sight passes through regions of enhanced 
density. Although the gaussian line shapes give a good fit 
to the spectra, we do not know the detailed velocity distri- 
bution of the shells along the line of sight, so we cannot re- 
liably determine the gas density contrast between the shell 
and inter-shell regions. However, the stronger components 
are clearly separate, and suggest that the contrast is a least 
a factor of a few. This is consistent with a large shell inter- 
shell contrast in the dust density derived from images in 
dust-scattered light by Mauron & Huggins (2000). 



Table 2. Parameters of the Ki components in IRC+10216 



Component 


Strength 
(rel.) 


b 

(kins" 1 ) 


Vr 

(kins" 1 ) 


1 


0.27 


2.8 


-32.4 


2 


0.32 


2.9 


-28.0 


3 


1.00 


5.5 


-21.0 


4 


0.51 


3.9 


-09.2 



3.2. Nai 

Fig. 4 shows the spectral regions covering the Nai D lines. 
The 5890 A line is in the upper panel and the 5896 A line 
is in the lower panel. The photospheric Na I lines of Star 6 
(marked with arrows at the right in each panel) are strong, 
and there is another weaker photospheric line in the lower 
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Fig. 4. Spectra of the Nai doublet at 5889.95 A {up- 
per panel) and 5895.92 A (lower panel). Arrows mark the 
photospheric and envelope components, as in Fig. 2. The 
triangles mark an interstellar Nai component (at V r ~ 
+9.5 kms- 1 ). 



panel. The dominant feature in the middle of each spec- 
trum is Nai absorption from the envelope of IRC+10216, 
centered near the systemic velocity. It can be seen that the 
circumstellar lines are extremely strong. All the compo- 
nents are highly saturated and the residual intensities are 
close to zero across both lines. 

There is an additional, weak Nai component that ap- 
pears in both spectra near +9.5 kms -1 (marked with trian- 
gles) which we identify as an interstellar component along 
this line of sight. Similar weak Nai absorption near this 
velocity is reported by Kendall ct al. (2002) along lines of 
sight to two other stars in this region of the sky, at angular 
distances of 153" and 2.5° from IRC+10216. ' 

We also searched for the much weaker ground state Na I 
line at 3302 A, which has an /-value ~ 70 times less than 
the 5890 A D line. The 3302 A line lies in a crowded, low 
signal-to-noise part of the spectrum, but the template ap- 
proach (Sect. 2) reveals a ~5 sigma detection, and the 
equivalent width is given in Table 1. The signal-to- noise 
ratio is too low to show any details of the profile. 

3.3. Ca I and Ca II 

Fig. 5 shows the spectral region around the Cai line at 
4226 A. The line falls in a very crowded region of the 
spectrum and provides a good illustration of the template 
method. 
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Fig. 5. Spectrum of the Ca I line at 4226.73 A. Upper panel: 
The observed spectrum (full line) and fitted template (dot- 
ted line). Arrows indicate the photospheric and envelope 
components as in Fig. 2. Lower panel: Spectrum normal- 
ized to the template. The dotted line replaces the spectrum 
over the photosperic core region. 



fa 




3967 



3968 



3969 



Wavelength (A) 

Fig. 6. Spectra of the Can H and K doublet at 3933.66 A 
(upper panel) and 3968.46 A (lower panel). Arrows mark 
the photospheric and envelope components, as in Fig. 2. 
The solid lines show the observed spectra, the dotted lines 
show the template spectra. 



The solid line in the upper panel shows the observed 
spectrum, and the envelope absorption is not immediately 
apparent. The dotted line in the upper panel shows the 
stellar template, scaled to match the observed spectrum. It 
can be seen that it gives a good overall fit to the data, and 
reveals the excess absorption from the envelope. The lower 
panel of Fig. 5 shows the envelope line, using the fitted 
template as the effective continuum. The small wavelength 
range around the core of the photospheric Cal line has 
been masked (with the horizontal dotted line) because the 
effective continuum is poorly defined at the low intensity 
of the line core. The main deviation from a flat continuum 
in the normalized spectrum is caused by slight differences 
between the Ca I line profiles of Star 6 and the Sun, which 
probably arise from slightly different surface gravities. 

It can be seen that the template fitting technique re- 
covers the Cal absorption in the envelope very effectively. 
The line shows a profile with narrow components similar to 
those seen in the K I lines. 

We also detected Can absorption in the H and K lines 
shown in Fig. 6. The 3933 A line is in the upper panel, 
and the 3968 A line is in the lower panel. For these spectral 
regions the crowding by photospheric lines is not severe, and 
the envelope absorption can be seen in the direct spectra. 
The adjacent Call photospheric lines are very strong, and 
their blue wings form the local (tilted) continuum for the 
circumstellar absorption. In the photospheric line core, the 



template is not a good fit on account of differences in the 
core reversals between Star 6 and the Sun. 

The Ca n line profiles show approximately the same pro- 
files as the Kl lines but are somewhat broader, by ~ 5- 
10 kms _1 (FWHM), suggesting additional absorption. This 
could be an interstellar contribution, or additional circum- 
stellar absorption. Compared to the neutral lines, the ion- 
ized Ca II lines sample the line of sight through the envelope 
to larger distances from the central star, and the kinematics 
of these outer regions have never previously been observed. 

3.4. Fei 

Figs. 7 and 8 show the spectral regions around the Fei 
lines at 3720 A and 3860 A, respectively. These fall in rel- 
atively crowded regions of the photospheric spectrum, and 
the template is needed to determine the envelope absorp- 
tion. The upper panel in each figure shows the observed 
spectrum (solid line) and the fitted stellar template (dot- 
ted line). The lower panel shows the normalized spectrum 
with the envelope absorption centered near the systemic 
velocity of IRC+10216. The signal-to-noise ratio of these 
spectra is lower than those discussed above, and the tem- 
plate fit is affected by differences in the solar and stellar 
spectra. The limited quality of the fit is probably responsi- 
ble for the fact that the equivalent width of the 3720 A line 
is slightly less than that of the 3860 A line, even though the 
f- value is larger (see Table 1). Nevertheless the lines are well 
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Fig. 7. Spectrum of the Fei line at 3719.93 A. Details as 
in Fig. 5. 



o 




3859 3860 
Wavelength (A) 

Fig. 8. Spectrum of the Fei line at 3859.91 A. Details as 
in Fig. 5. 



detected, and the multi-component character of the enve- 
lope absorption is similar to that seen in the lines of the 
other species. 

We also searched for two additional Fe I lines at 3824 A 
and 3441 A. The 3824 A line is weaker on account of a sig- 
nificantly lower /-value, and is marginally detected at ^4 
sigma. The 3441 A line lies at shorter wavelengths where 
the spectrum is poorer. It is detected with a strength com- 
parable to the Fei lines shown in Figs. 7 and 8, but with a 
lower signal-to-noise ratio. 

3.5. Other lines 

Searches were also made for lines of several other species, 
including Al I, Mn I, Cr I, Ti I and Ti II, and Sr II, as listed in 
Table 1. With the exception of Cn, no significant envelope 
absorption was detected. For Cn, weak (^2-3cr) features 
are seen at the wavelengths of most of the accessible lines. 
When appropriately weighted by the oscillator strengths 
and noise levels they provide an overall 4-sigma detection 
of this species. 



4. Abundances in the envelope 

4.1. Column densities 

The equivalent widths (W\) of the absorption lines ob- 
served in the envelope are listed in column 4 of Table 1. 
For the non-detected lines, we list 3-sigma limits given by: 



W x < 3crAA/V^, 



where a is the noise level and n is the number of resolution 
elements in the line width AA. 

The column density (N) derived for each species is given 
in column 5 of Table 1. It is based on the equivalent widths 
and the /-values given in the table. For the weak lines, the 
column density is directly related to the equivalent width 
by the optically thin formula: 

N = 1.13 x 10 20 Wx/fX 2 , 

where Wx and A are in A units. 

For the stronger lines, where the optical depths are 
larger, we determine the relation between Wx and N using 
the multi-component model used to fit the high signal-to- 
noise profiles of the Ki lines (Table 2). With this model 
we find that the Ki lines are moderately saturated; the 
column densities derived from the 7699 and 7665 A lines 
are larger than the optically thin values by factors of 1.7 
and 2.3, respectively. The column densities from the two 
lines agree to within 20%, and are consistent with a limit 
of iV(K i) ^ 3x 10 13 cm~ 2 determined from the upper limit 
to the optically thin 4044 A Ki line. For Nai, the D-lines 
are extremely optically thick, and therefore insensitive to 
the column density, although the very low residual inten- 
sity across the lines yields a lower limit of iV(Nal) <;1- 
2 x 10 13 cm~ 2 . Fortunately, we also detect the optically 
thinner Nai line at 3302 A line, and we use this with the 
multi-component saturation curve to determine the column 
density of Nai given in Table 1. The result is consistent 
with the lower limit from the D-lines. 

The profiles of the ionized lines might be expected to 
differ from those of the neutral lines because the absorption 
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can occur at different locations along the line of sight. This 
does not affect the column density estimates for Srii and 
Tin because the lines are very weak (below the detection 
level) and are therefore optically thin. For the stronger Ca II 
lines, the profiles are found to be similar to the K I profiles, 
but as noted in Sect. 3.3, they are slightly broader with 
an additional contribution to the absorption. We therefore 
model the Call lines as the sum of two contributions; the 
fit of the K I multi-component model to the main part of 
the profile, and an additional, optically thin component to 
give the total equivalent width. The additional component 
contributes 32% and 17% of the equivalent widths (22% 
and 13% of the derived column densities) for the 3933 and 
3968 A lines, respectively. We assume that the additional 
absorption is circumstellar, but since it is relatively small, 
the results do not depend sensitively on this assumption. 

4.2. Column density of hydrogen 

Although there are no direct observations of hydrogen along 
the line of sight, we can use estimates of the mass loss rate 
of the envelope to determine the hydrogen column density. 
Since the early work by Kwan & Hill (1977) there have 
been numerous estimates of the mass loss rate based on 
millimeter CO observations, and the results are fairly con- 
sistent when the different distances, CO abundances, dust- 
gas heating rates, and He content are taken into account. 
For a distance of 120 pc, we adopt a mass loss rate Mr 
(in hydrogen) of 1.25 x 10~ 5 M Q yr _1 (corresponding to a 
total mass loss rate of 1.75 x 10~ 5 M yr -1 ), based on the 
analysis of Schoier & Olofsson (2001) approximately cor- 
rected for the effects of He. This value is consistent with 
other recent estimates. 

For A#h = 1-25 x 10 -5 M Q yr _1 , the column density 
of hydrogen N(R), where iV(H) = TV (Hi) + 27V (H 2 ), is 
1.3 x 10 21 cm~ 2 along the line of sight 35" from the cen- 
ter. The uncertainty is a factor ~ 2, which results from 
uncertainties in M. In addition to a smooth decrease of 
the column density with distance from the center, there are 
other variations caused by the multiple shell structure in 
the envelope. From an analysis of scattered light images 
(e.g., Fig. 1) we find typical variations in the dust column 
density of ±20%, and somewhat smaller variations near the 
region of Star 6. These variations are much smaller than 
the density contrast of the shells because the column den- 
sity averages the density along the line of sight through 
several shell and inter shell regions. The envelope structure 
is therefore not a major source of uncertainty for N(H). 

4.3. Ionization fraction 

The line of sight passes through regions of the circumstellar 
envelope where the metals are partially ionized. Since most 
of the observed lines arise from single stages of ionization, 
we need to consider ionization corrections in order to deter- 
mine the total gas phase abundances. We estimate the cor- 
rections from the relative column densities of the ionization 
stages of each metal, obtained using the photo-ionization 
model of an expanding envelope discussed by Glassgold & 
Huggins (1986). 

In the model, neutral atoms (which may result from the 
dissociation of molecules) emerge from the dense, shielded, 
inner envelope and are photo-ionized by the ambient inter- 



stellar radiation field. We calculate the ionization fraction 
as a function of radius, and use this to calculate the ioniza- 
tion fraction along the line of sight. The ionization of each 
element is governed by equation 4.9 of Glassgold & Huggins 
(1986), with no contribution from chromosphcric radiation. 
The ionization depends on the photo-ionization rate (given 
by the interstellar rate and the envelope shielding) and re- 
combination. In solving the ionization equations we use the 
interstellar photo-ionization rates and recombination rates 
from Pcquignot & Aldrovandi (1986), except for Cr and Sr 
which are not included in their compilation; for these the 
photo-ionization rates are from Glassgold & Huggins (1986) 
and the recombination rates from Bernat (1976). The ra- 
dial dependence of the electron abundance is adopted from 
Cordiner et al. (2007) , although recombination is relatively 
unimportant except for Al. For the shielding of interstellar 
radiation in the envelope we use the standard carbon dust 
model and dust-to-gas ratio of Cherchneff ct al. (1993) and 
the shielding function of Morris & Jura (1983). 

There are a number of uncertainties in the parameters 
of the model, the most important being the strength of the 
ambient radiation field (which determines the interstellar 
photo-ionization rates) and the shielding in the envelope 
(which includes uncertainties in the dust parameters, the 
mass loss rate, and the envelope geometry) . Fortunately we 
observe both the neutral and ionized column densities of 
Ca, and we use this to set the ionization level in the model. 
For the nominal parameters given above, the predicted ra- 
tio iV(Caii)/A(Cai) is a factor 3.6 larger than observed. 
This is fair agreement considering that the ionization in 
the envelope has not previously been constrained in this 
way. However, we can do much better, by fine tuning the 
ionization level in the envelope to fit the Ca ionization ex- 
actly. This can be done by adjusting the ambient radiation 
field (by a factor of 0.46) or the shielding optical depth (by 
a factor 1.5). Both variations are within their respective 
uncertainties. For specificity, we adopt the reduced radia- 
tion field to calculate the ionization of the other elements. 
The resulting ionization corrections (Ci = A^(i + n)/iV(l) 
for neutral species, and A^(i + n) /N(u) for singly ionized 
species), are given in column 4 of Table 3. Adjusting the 
shielding instead, by the amount given above, produces es- 
sentially the same ionization corrections. 

The ionization correction for All in Table 3 is much 
larger than for the other metals because of its relatively 
high photo-ionization rate. Our observations are therefore 
not very sensitive probes of the column density of Al in 
the gas phase because the Al is nearly completely ionized 
along the line of sight. The ionization corrections for the 
other neutral species are much smaller. For example, Na 
is predominantly neutral. This is in contrast to typical in- 
terstellar clouds with a similar column density of hydro- 
gen. In the circumstellar envelope, the characteristic out- 
flow time at 35" is shorter than the photo-ionization time 
given by Glassgold & Huggins (1986). Hence, even with- 
out dust shielding, the Na atoms are expected to be largely 
neutral, as given in the table. 

4.4. Abundances 

Except for Ca and Ti, the total gas phase column densities 
of the metals are determined from the observed column den- 
sities and the ionization corrections. The results are given in 
column 5 of Table 3. Although the iV(Caii)/./V(Cai) ratio 
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was used to determine the ionization level in the envelope, 
we use the sum of the ionization stages to obtain the total 
Ca column density, independent of the ionization. Similarly, 
the limit for Ti is based on the observed limits for Ti I and 
Tin, and so is independent of the ionization. 

The gas phase abundances of the metals relative to hy- 
drogen (X) are determined from the total column densities 
and the value of iV(H) from Sect. 4.2, and are given in col- 
umn 6 of Table 3. For reference, the solar values (Xq) are 
given in column 7 of the table, taken from Lodders ([2003). 
Comparison of the envelope and solar abundances shows 
that there are large deficiencies in the gas phase abundances 
in the envelope, which vary from metal to metal. 

5. Dust condensation 

Several processes affect the state of the circumstellar ma- 
terial as it moves from the stellar photosphere out through 
the envelope into the interstellar medium. In order of in- 
creasing distance from the star these processes include: 
dust condensation in conditions of approximate thermo- 
dynamic equilibrium; gas phase chemical reactions; photo- 
dissociation of the molecules; and eventual photo-ionization 
of the atomic constituents in the outer envelope (e.g., 
Gilman 1969; Tsuji 1973; McCabe et al. 1979; Huggins & 
Glassgold 1982; Lafont et al. 1982). In the carbon-rich en- 
vironment of IRC+10216, the main component of the dust 
is amorphous carbon with a minor component of SiC (e.g., 
Martin & Rogers 1987). The state of other elements, es- 
pecially the refractory metals, is not well understood (e.g., 
Turner 1995). Hence our observations of the gas phase met- 
als in the outer envelope provide new constraints on the gas 
phase chemistry and the dust condensation. 

5.1. Observed depletions 

Comparison of the envelope abundances with the solar 
abundances in Table 3 shows that most of the gas phase 
metal atoms in the envelope are "missing" . It is most un- 
likely that they are in the form of gas phase molecules. The 
fractional abundance of even the most extreme case of Ca, 
where X(Ca)/ Xq ~ 3 x 10 -3 , is an order of magnitude 
larger than the largest fractional abundance of any metal 
bearing molecule detected outside of the core region (see 
Table 4). In addition, most molecules are dissociated closer 
to the star than the 35" offset of our line of sight. It is 
therefore reasonable to infer that the atomic abundances 
observed are good approximations to the total gas phase 
metal abundances in the envelope, and that the missing 
atoms are depleted onto dust grains. 

In column (8) of Table 3 we give the conventional mea- 
sure of depletion log 10 8, where 8 = X/Xq is the depletion 
factor. The observational limits for the abundances of Al 
and Sr are <^ the solar values, so in these cases we have no 
significant limits for <5, although Sr is an s-process element 
and may be enhanced in IRC+10216 and depleted. 

Based on the measured depletions, there are two im- 
mediate conclusions. First, the metals in this carbon-rich 
archetype are primarily in the form of solids, and this is 
the dominant form returned to the interstellar medium. 
Second, in spite of the depletion, a significant residue of 
metallic atoms remains in the gas phase and varies from 
metal to metal. 



5.2. Condensation and adsorption 

There are no detailed predictions for the depletion of met- 
als in IRC+10216, but there are some important consider- 
ations that bear on the issue. A commonly used approach 
to the condensation of solids in circumstellar envelopes is 
the assumption of thermodynamic equilibrium in the warm, 
dense, inner envelope, where the chemical time scales are 
rapid compared with the expansion time scale. Under these 
conditions the formation of solid particles is controlled by 
the condensation temperature of the primary condensate of 
each species. 

For a carbon-rich envelope, the condensation sequence 
depends somewhat on the C/O ratio and the gas pressure. 
The following sequence, for C/O = 1.1 and a pressure of 
10~ 6 bar (from Lodders & Fegley 1995, updated for Fe by 
Lodders & Fegley 1999) is representative: C (1670 K), TiC 
(1640 K), SiC (1460 K), FeSi (1230 K), A1N (1170 K), CaS 
(1150 K), MgS (960 K), with other more volatile metals 
such as Na and K at lower temperatures. Cr and Mn prob- 
ably form sulphides but their location in the sequence is 
uncertain. Thus the metals Ti, Fe, Al, Ca, and then K and 
Na, are expected to be removed from the gas phase succes- 
sively. Those with lower condensation temperatures are less 
likely to go to completion because of the decreasing density 
with radius. 

This qualitative picture is largely consistent with the ob- 
served depletion pattern. The observations provide a firm 
upper limit on the gas phase abundance of Ti (from the Ti I 
and Till lines); Fe and Ca are strongly depleted, although 
Ca is more depleted than Fe (in reverse order to the con- 
densation sequence); and Na and K are less depleted. 

The simplifying assumption of thermodynamic equilib- 
rium is not a complete physical theory of condensation be- 
cause kinetic effects must play a role near freeze-out. In 
addition, once formed, the grains can act as sites for the 
adsorption of gas phase species further out in the enve- 
lope (Jura & Morris 1985). The adsorption depends on the 
binding energy of the species to the grain surface, and on 
the sticking probability (p), which is essentially unknown. 
Using the analysis of Jura & Morris (1985), Turner (1995) 
finds that for IRC+10216, a volatile species such as K is 
depleted by adsorption to a fractional abundance of 0.84 
(for p — 0.1) and 0.17 (for p = 1), and a more refractory 
species such as Al is depleted to a fractional abundance of 
0.05 (for p = 0.1) and lO" 10 (for p = 1). Thus adsorption 
of metals may be as important as the initial condensation. 

Our observations of a residual atomic component in the 
gas phase constrain the efficiency of both condensation and 
adsorption. Below about 1150 K, the phase diagrams of 
Lodders & Fegley (1995) show that Ca is essentially com- 
pletely removed from the gas phase, but our observations 
show that even for this refractory metal there is a residual 
component in the gas phase. Metal bearing molecules with 
Na and K (albeit with low fractional abundances) are seen 
in the core region (see Sect. 6) but not in the extended en- 
velope. This may be the result of adsorption. On the other 
hand, adsorption of the more refractory species to levels of 
10 -10 are clearly ruled out by the observations. 

We conclude that current theoretical ideas are qualita- 
tively in accord with our observations of metal depletion, 
but improved models with specific quantitative predictions 
are needed to discriminate the underlying processes. 
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Na 


I 


4.6(14) 


1.22 


5.6(14) 


4.2(-7) 


2.00(— 6) 


—0.68 


—0.06 


—0.95 


A 1 

Al 


1 


<2.8(12J 


3.5(3) 


<9.8(15) 


<7.3(— 6) 


2.88(— 6) 
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K 


I 


4.6(12) 


1.60 


7.4(12) 


5.5(— 9) 


1.29(— 7) 


— 1.37 


—0.17 


—1.09 


Ca 


1 


1.9(12) 


4.70 


8.9(12) b 


6.6(-9) 


2.19(-6) 


-2.52 


-0.75 


-3.73 




II 


7.0(12) 


1.27 














Ti 


I 


<1.9(12) 


4.26 


<3.9(12) b 


<2.9(-9) 


8.32(-8) 


<-1.45 




-3.02 




II 


<2.0(12) 


1.31 














Cr 


I 


1.4(12) 


10.6 


1.5(13) 


l.l(-8) 


4.47(-7) 


-1.60 




-2.28 


Mn 


I 


<3.1(12) 


2.21 


<6.9(12) 


<5.1(-9) 


3.16(-7) 


<-1.79 




-1.45 


Fe 


I 


8.8(13) 


2.58 


2.3(14) 


1.7(-7) 


2.95(-5) 


-2.24 


-1.68 


-2.27 


Sr 


II 


<2.7(11) 


3.82 


<1.0(12) 


<7.6(-10) 


8.13(-10) 


C 







a Ionization correction d = N(I + II)/N(I) for neutral species and N(I + II)/N(II) for singly ionized species. 
b Based on observed values of N(l) + N(H), independent of d. 
c Upper limit on X > Xq, no useful limit on S. 



5.3. Comparison with PNe 

Carbon-rich planetary nebulae (PNe) are the immediate 
descendants of carbon-rich AGB stars, in which the cir- 
cumstellar material has undergone major changes. The gas 
has been photo-ionized, and the dust grains are exposed to 
intense radiation fields and high (~ 10 4 K) temperatures. 
A comparison of depletions in AGB stars and PNe may 
therefore reveal some aspects of grain evolution. 

Element abundances have been extensively measured 
in the ionized gas in PNe. The abundances are, however, 
subject to systematic, and sometimes, large uncertainties, 
and relatively few metals have accessible lines covering 
the appropriate stages of ionization. For comparison with 
IRC+10216 we focus on NGC 7027, which is one of the 
most intensively studied, carbon-rich PNe. NGC 7027 is 
relatively young and still surrounded by a substantial enve- 
lope of molecular gas (Cox et al. 2002 i. The circumstellar 
conditions before the formation of the nebula were therefore 
somewhat similar to the current state of IRC+10216. 

There have been numerous abundance analyses of 
NGC 7027. We have taken the abundances of Na, K, Ca, 
and Fe from recent, comprehensive studies by Keyes et 
al. (1990), Middlemass et al. (1990), Bernard Salas et al. 
(2001), and Zhang et al. (2005), and we give the correspond- 
ing depletions in column (9) of Table 3. For species in com- 
mon the depletions have been averaged. Even with these 
state-of-the-art analyses, the differences in abundances be- 
tween the different studies range up to a factor of ~ 3. 

Comparison of the depletions in IRC+10216 and 
NGC 7027 reveals some significant differences. First, the 
metals Na and K are much less depleted in the ionized 
nebula, where they are close to the solar values. Second, Fe 
and Ca are still significantly depleted in the nebula but, evi- 
dently less than in the circumstellar envelope. These results 
suggest the following evolutionary effects in the transition 
from AGB star to PN: the nearly complete evaporation of 
the volatile species Na and K, and the partial erosion of 
more refractory species Ca, and possibly Fe. Further evi- 
dence for this view comes from the study of NGC 7027 by 



Kingdon et al. (1995), who argue that Ca is depleted by 
more than 2 orders of magnitude near the periphery of the 
nebula (as in IRC+10216), but is much less depleted near 
the center of the ionized nebula. 

The abundances in NGC 7027 are fairly similar to other 
carbon-rich PNe, e.g., the Fe depletion is similar to the 
typical Fe depiction found in a sample of low ionization PNe 
by Delgado Inglada et al. (2009). The trends noted here 
may therefore be a general characteristic of the evolution 
of dust from the AGB to PNe. 



5.4. Comparison with the ISM 

It is also of interest to compare our results for IRC+10216 
with depletions in the ISM. We take the line of sight towards 
£ Oph as representative of the ISM, bearing in mind that 
the overall level of IS depletions varies with location but the 
general pattern remains the same. The depletions towards 
£ Oph, from Savage & Sembach (1996), are listed in column 
(10) of Table 3. 

In spite of the different physical and chemical environ- 
ments that lead to grain formation in IRC+10216 and the 
ISM, it can be seen that the depletion patterns are qualita- 
tively similar. Na and K are the least depleted, at compara- 
ble levels; Mn, Fe, and Cr are probably similar although the 
detailed pattern may differ; and Ca is the most depleted in 
both data sets. 

In studies of the ISM it is found that element depletions 
correlate with the (oxygen-rich) condensation temperature 
of the clement-bearing solid. This was interpreted in terms 
of grain formation at high temperatures in the winds of 
mass- losing giants (Field 1974). More recent studies indi- 
cate that a significant component of the dust in the ISM is 
formed in situ. The correlation with the condensation tem- 
perature may therefore reflect some aspect which is shared 
by the formation process in the ISM. The similarity of the 
pattern that we find in IRC+10216 may have some bearing 
on this question, and deserves further attention. 
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Table 4. Metal bearing molecules in IRC+10216 
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NaCl 


2.5 


1.3(-09) 


6.3(-4) 


1 


NaCN 


2.5 


6.8(-09) 


3.4(-3) 


1 


KC1 


2.5 


3.8(-lQ) 


2.9(-3) 
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A1F 


2.5 


4.6(-08) 


1.6(-2) 


1 


A1C1 


2.5 


2.8(-08) 


9.8(-3) 


1 


A1NC 


5-15 


2.3(-10) 


7.9(-5) 


3 


MgNC 


10-20 


6.3(-09) 


1.8(-4) 


1 


MgCN 


10-20 


2.9(-10) 


8.0 (-6) 


4 



a Source radius. 

b Abundance relative to H, for d = 120 pc, M(H) = 1.25 x 
10" 5 M Q yr-\ 

References: (1) Highberger & Ziurys (2003). (2) Assuming 
X(KC1) = 0.3X( NaCl) Cernicharo & Guelin HT557> . 
(3) Ziurys et al. 1(1002)1 . ( 4) A ssuming A(MgCN) = 
0.045X(MgNC) Ziurys et al. (JTJ5S]). 



6. Metal chemistry 

The circumstellar envelope of IRC+10216 exhibits a re- 
markable gas phase chemistry, with more than 50 molecular 
species detected to date, mainly through their rotational 
transitions at millimeter wavelengths. The majority of the 
molecules are formed of the abundant elements H, C, N, 
and O, several include Si and S, but a small number un- 
expectedly include metals, Na, K, Mg, and Al. The first of 
these were discovered by Cernicharo & Guelin (1987). Our 
observations of gas phase metal atoms in IRC+10216 con- 
strain some aspects of the metal chemistry in the envelope. 

The metal bearing molecules detected in the envelope 
are listed in Table 4. In addition to these, numerous other 
metal bearing species have been searched for at comparable 
levels, but not detected (e.g., Turner 1995). There are also 
other metal species whose rotational frequencies have not 
yet been measured in the laboratory (e.g., Ziurys 2006b). 

For comparison with the atomic abundances reported 
here, we list updated abundances for the metal molecules 
in column (3) of Table 4. These are based on the column 
densities given in the references in the table, but are de- 
rived for the distance and mass loss rate used in this pa- 
per, and are relative to hydrogen. The fractional abun- 
dance X mo \/ Xq, is the fraction of each metal in a particular 
molecular species. 9 r is the measured, or inferred, angular 
radius of the molecular distribution, adopted from the ref- 
erences. 

The metal bearing molecules divide naturally into two 
groups. The first is confined to the dense core region repre- 
sented by 6 T ^ 2.5" in the table. The second, which includes 
only Mg and Al bearing molecules, is found in distinct shells 
in the envelope, in the photo-dissociation region. The gen- 
eral level of incorporation of the metals into molecules is 
very small in both groups. 

In the outer envelope, the fractional abundances of the 
metal bearing molecules are ^ 0.02 per cent, while the 
atoms that we detect have fractional abundances of 0.3- 
20 per cent. This we used in Sect. 5.1 to justify the as- 
sumption that the observed atoms are the dominant gas 
phase metal carriers along the line of sight. 



Even in the core region the molecules do not seem to 
be dominant. The observed fraction of Na in the core in 
the form of molecules is ~ 0.4 per cent, and we detect 
20 percent in the form of atoms farther out. Similarly, the 
fraction of K in the core in the form of molecules is ~ 0.3 per 
cent, and we detect 4 per cent in atoms farther out. For 
these species the numbers are consistent with the idea that 
apart from condensation or adsorption onto dust grains, 
the dominant gas phase species is atomic throughout the 
envelope. There is lack of specific molecular information 
for the other metals, but the absence of detected molecular 
species argues that this applies to them as well. 

Our finding that residual metals are present in the enve- 
lope in the form of neutral atoms and ions provides an ob- 
servational foundation for understanding certain aspects of 
the metal chemistry. For example, it has been proposed that 
metal cyanides are formed by reactions of metal ions with 
cyanopolyyenes (Dunbar & Petrie 2002), and the model 
of Cordiner & Millar (2009) shows that this can account 
for the observed abundance of MgCN if enough Mg + is 
present in the gas phase. We have not observed Mg, but 
we expect that it behaves like the other metals. Its con- 
densation temperature is significantly less than that of Ca 
and Fe, so it is likely to be less depleted in the envelope. 
The gas phase abundances of the metals that we observe 
are given in column (6) of Table 3. These may undergo re- 
actions with cyanopolyyenes and other abundant neutral 
molecules, such as unsaturated hydrocarbons, and lead to 
a variety of metal bearing species. The largest abundances 
that we observe in the gas are those of Na and Fe (whose 
large depletion is balanced by its high cosmic abundance). 
Thus Fe offers interesting possibilities for a potentially ob- 
servable Fe chemistry. 



7. Conclusions 

The observations of IRC+10216 reported in this paper rep- 
resent the first comprehensive study of atomic metals in a 
carbon-rich circumstellar envelope. We detect lines of Na, 
K, Ca, Cr, and Fe, and obtain upper limits for Al, Ti, Mn, 
and Sr. Combined with a simple model of the ionization, 
the observations provide estimates of the gas phase metal 
abundances in the outer envelope. 

The results show that the metals, especially Ca and Fe, 
are significantly depleted onto dust grains in the circum- 
stellar envelope, and this is the dominant form returned to 
the ISM. The depletion pattern has some similarity with 
depletion in the ISM, and is roughly consistent with expec- 
tations of dust condensation in a carbon-rich envelope. 

Although the metals are depleted in the envelope, 
atomic metals in the form of neutral atoms and ions ap- 
pear to be the major metal species in the gas phase. As 
such, they likely play a key role in the metal chemistry of 
the envelope. 
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